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SUMMARY 


The stability of an electrically conducting fluid subjected to two 
dimensional disturbance was investigated. The physical system consists of two 
parallel infinite vertical plates which are thermally insulated. Applied 
normal to the plates is an external magnetic field of constant strength. The 
fluid is. heated from below so that a steady temperature gradient is maintained 
in the fluid. The governing equations were derived by perturbation technique, 
and solutions were obtained by a modified Galerkln method. It was found that 
the presence of the magnetic field increases the stability of the physical 
system, and instability can occur in the form of neutral or oscillatory 
instability. 


INTRODUCTION 


The Bernard problem of stability of a viscous, stationary fluid heated 
from below has been a subject of many investigations (ref. 1). It is found 
that when the temperature gradient exceeds a critical value, instability sets 
in as stationary cellular convection. For the case of an electrically con- 
ducting fluid in the presence of a magnetic field, however, the physical 
phenomenon becomes more complicated. The motion of the fluid crossing the 
magnetic field, causes electrical currents to be generated, and the current 
carrying fluid elements traversing the magnetic field gives rise to an 
additional body force (the Lorentz force) which tends to retard the fluid 
motion. Consequently, the stability of the system is greatly increased. 
Furthermore, Chandrasekhar (ref. 1) proved that depending on the relationship 
between the electric resistivity and the thermal diffusivity of the fluid, 
instability can manifest itself in the form of stationary or oscillatory 
motion. 

Exact solution of instability due to small two dimensional disturbance of 
a viscous fluid bounded by two parallel vertical planes was obtained by Ylh 
(ref. 2). The bounding planes are thermally insulated, and ah upward tempera- 
ture gradient is applied to the fluid. Ylh showed that for disturbances 


1509 



periodic in the vertical direction, the most unstable modes are associated 
with the wave number zero. It is also shown that instabilities due to anti- 
symmetric disturbances with respect to a median vertical plane are more easily 
excited than those due to symmetric disturbances of the same wave number. 

In this investigation, it is proposed to extend Ylh's finding by 
incorporating the effect of a magnetic field applied normal to the bounding 
walls. The governing equations were simplified by using Boussinesq approxi- 
mation. The onset of neutral as well as oscillatory instabilities was 
studied. 


SYMBOLS 

d distance between the vertical planes 

g acceleration due to gravity 

intensity of applied magnetic field of constant strength 
h perturbation quantity of magnetic intensity 

n wave number 

P pressure 

Pr^ Prandtl niomber, v/k 

Pr 2 magnetic Prandtl number, v/n 

f 

Pr^ parameter defined as p/K 

2 2 

Q electromagnetic number, (pH d )/(4irpvp) 

, 2 2 2 
Q parameter defined as (pH d )/(4irpK ) 

R Rayleigh number, (ga3d^)/(Kv) 

u,w perturbation velocity quantity in x and z direction 
x,z space co-ordinates 

a coefficient denoting change of fluid density per degree rise in 

temperature 

g temperature gradients in fluid 

n electric resistivity 

9 fluid temperature 
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y 

V 

p 

0) 

c 

6o3 


thermal diffuslvlty 
magnetic permeability 
kinematic viscosity 
density 

critical wave speed 

parameter defined as (6P/p) + (UH^h^) / (4 ttp) 


Laplace operator, — ^ 

3x 


3z' 


Superscript: 

* dimensionless 

Subscripts: 

x,z X and z component, respectively 
c critical 


FOiJMULATION 


The physical system under study consists of two infinite and thermally 
insulated vertical planes placed at a distance d apart. Applied normal to the 
bounding planes is an external magnetic field of constant strength Hq. The 
electrically conducting fluid is heated from below so that a steady tempera- 
ture gradient is maintained in the fluid. Figure 1 shows the schematic 
diagram of the physical system and the Cartesian co-ordinate adopted for this 
study. 


A detailed development of the governing differential equations and the 
perturbation equations are presented in chapters 2 and 4 in reference 1. For 
the present problem, the perturbation equations are given below: 

Continuity equation 


3u ^ 

9x 3z 


0 


Momentum equation 


iH. 

9t 


3x 


- ^®o 

((Su)) + vV^u + 



3x 


( 1 ) 


( 2 ) 
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3w 

3t 

MH 3h 

= - (6a,) + gae + vV^w + 3^ 

(3) 

Energy equation 



30 

3t 

= gw + kV20 

(4) 

Equations for E-M 

field 


3h 

X 

3x 

3h 

+ ^ = 0 
3z 

(5>^ 

3h 

X 

3t 

= H 1^ + nv^h 
0 3x X 

(6) 

3h 

z 

3t 

= E ^ + nV^h 
o 3x z 

(7) 


The physical variables used In the above equations are defined In the SYMBOLS, 
In order to reduce the above equation to dimensionless form, the following 
dimensionless quantities are Introduced: 


* ud 
u = — 


^ = d ’ 


ft* = ^ 
gd 


* wd 

w = — > 


* X 
X = 


d ’ 

. h 4. h 

u* - ^ u* = _z 


H 


* 

t = 


o 

Kt 


H 


-* d — 

6(jo = — 2 <5u 

K 


By employing these dimensionless quantities, eqs. (1) to (7) can be 
expressed in non-dimension form as in the following: 

* 


3u 

■f 

9x 


3u 

3t 


3w 

* 

3t 


36 

3t 


+ 


3w 

■f 

3z 


= 0 


( 5 '^) + Pr V*^u* + q' 
3x 


3h, 


3x 


^ (6(jj ) + Pr^R 0 + Pr^^V ^w + Q 


3z 


* *2 * 
= w + V 0 



(8) 


(9) 

* 


3h 
z^ 

A 

(10) 

3x 

(11) 


9h 


X 


3h 


3z 


= 0 


( 12 ) 
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3h 


3t 


9u 

* 

3x 


* * 

+ Pr- V h 

2 X 


( 13 ) 


3h 


3t 


3w 

i 

3x 


* *2 * 
+ Pr_ V ^ h 
^ z 


(14) 


Elimination of 3o) from eqs. (9) and (10) yields 

[1— - Pr^^ V*2] (lu_ _ 

3t 3z 3x 3x 

, „ 3h* 3h* 

+ Q 

3x 3z 3x 


(15) 


The system of equations can be simplified by introducing the following 
stream functions ip and <|) which satisfy eqs. (8) and (12) automatically. 


u 


= - 


3z 

h * = - ^ 


8x 


(16) 


h * = M. 


X - z 

3z 3x 


Substituting eq. (16) into eqs. (15), (11) and (14), the governing 
eqs. become 

- Pr V*2j ^ q' ^ (V%) (17) 


3t 

.3 


'1 * 
3x 


3x 


[ * 
3t 


- e* . ^ 

9x 


__ Pr ’v*2] M_= 3JL_ 
* ^‘'2 •* * *2 
3t 3x 3x 


(18) 

(19) 


at 


Eqs. (17), (18), and (19) are subjected to the boundary conditions 
X* = + y. Ip = Dij) = D9* = <t> = 0 (20) 


SOLUTION 


Solutions were sought by using Galerkin method modified by Finlayson 
(ref. 3). Assuming disturbances of the following form 

ip = A(t) ip(x) e^*^^ 

e* = B(t) 0(x) e^“^ (22) 
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<|) = C(t) $(t) e 


inz 


(23) 


In selecting an appropriate trial function for a disturbance 
symmetrical with respect to the medium plane was chosen as 


^ _ Cosh X^X 

i 


Cos X^X 


Cosh (X^/2) Cos (X^/2) 
where the characteristic values X. are given by 


(24) 


tanh (X^/2) + tan (X^/2) = 0 
The assumed function for Q and s are 

* 

9 = sin a^x , s = sin 3^x 

where “i = (2j-l) tt , 3^ » 2jir and j = 1,2,3, etc. 


(25) 


(26) 


Substituting the trial functions given by eqs. (24) and (26) into the 
governing eqs. (17), (18) and (19), and after diagonalization, the governing 
differential equations are reduced to 


dA 


[(C,/C, - n"(C,/C,)] PrT[(C,/c/^) - 2n^(C,/C,. ) 


3 1 


dt 


3 1 


3 1 


/ ^ 

+ n (C./C.)] A + Pr^Ra. (C./Cosa.X )B + 

Ix j 1 

q' [-3^(6^^^ + n^)(Cj/Cos3^x*)]C (27) 

dB ^ 2 2 

(sina ,/sina ,) -j- = (sina,/C.)A+ [-(a +n ) (sina./sina.) ]B (28) 

3 X dt XX j X 

(sin3j/sin3.) (sin3 . /C . ) A + Pr„ [ (3 . ^+n^) (sin3 . /sin3 . ) ]C (29) 

j»xai: J X 

Where the inner products are defined in the following manner 

1 

;,dx* (30) 


(C 


,/c/') - c.o\.t 

’ JLx ^ " 


The method for determining the stability criteria for eqs. (28), (29) 
and (30) is presented in reference 3. 


DISCUSSION OF RESULTS 


Numerical results were obtained by using a digital computer. The critical 
Rayleigh number was obtained for different combinations of the Prandtl number 
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Pr^, the electromagnetic Prandtl number Pr 2 , the electromagnetic number Q and 
wave number n. The results were presented in graphical form from figures 2 to 

5- 


For the case when Q = 0, the present problem was reduced to the Ylh's 
problem. Examination of the results presented in Table 1 shows that the 
critical Rayleigh number converges rapidly to Ylh's exact solution with the 
second approximation, thus confirming the validity of the results. 

Table 1: COMPARISON OF R FOR Q = 0 

c 



n 

0 

1 

2 

4 


Yih's soln. 

500.8 

592.0 

856.0 

2952 

R 

c 

1st approx. 
2nd approx. 

514.6 

501.0 

597.1 

578.8 

885.5 

857.0 



3100 

2927 


For the case of neutral stability, the critical Rayleigh number is only a 
function of the wave number n and Q which is a measure of the strength of the 
applied magnetic field (see figs. 2 and 3). The Prandtl numbers Pr^ and Pr 2 
have no effect on the stability of the system. This is expected since for 
the case of neutral stability, the terms Involving the time derivative in the 
governing equations (eqs. 28 to 30) are dropped, and Pr^^ and Pr 2 would be 
cancelled from the resulting equations. Similar situations arise in the 
Bernard problem. 

For the present physical system, overstablllty can occur. Figure 4 shows 

the variation of R as a function of Q, Pr. and Pr.-,. The results show that an 

c 12 

increase in either Pr^ or Pr 2 tend to decrease R^. This finding is in agree- 
ment with the results obtained for Bernard Problem (Chapter 4, ref. 1). 

With the trial functions used in this study, the convergence of the 
numerical results are fairly satisfactory. This is especially so in the case 
of neutral stability; an average discrepancy of about 5% was observed between 
the first and second approximation. For the case of overstability, greater 
deviation is observed. A difference of 25% was observed between the first and 
second approximations. Higher order approximation would be needed to obtain 
more accurate results. 


CONCLUDING REMARKS 


The following conclusions can be drawn from this study: 

1. The presence of a magnetic field Increases the stability. 
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2. In the presence of a magnetic field. Instability can occur in the 
form of neutral or oscillatory instability. When oscillatory instability 
occurs, the effect of Pr^^ and becomes important. An Increase in either 
one or both parameters tend to render the system less stable. 

3. The most unstable modes are associated with the wave number zero. 
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U) for onset of overstablllty . 



